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Abstract 

A quantum system S interacts in a successive way with elements f of a chain of 
identical independent quantum subsystems. Each interaction lasts for a duration 
T and is governed by a fixed coupling between S and £. We show that the system, 
initially in any state close to a reference state, approaches a repeated interaction 
asymptotic state in the limit of large times. This state is r-periodic in time and 
does not depend on the initial state. If the reference state is chosen so that S 
and £ are individually in equilibrium at positive temperatures, then the repeated 
interaction asymptotic state satisfies an average second law of thermodynamics. 



1 Introduction 

In this introduction we outline our main results and the relevant ideas of their proofs. 

Suppose a quantum system S interacts with another one, £, during a time interval 
[0, r), where r > is fixed. Then, for times [r, 2r), S interacts in the same fashion with 
another copy of £, and so on. The assembly of the systems £ (which are not directly 
coupled among each other) is called a chain, C = £+£ + ■■■ . The system S + C, with an 
interaction as described above, is called a repeated interaction quantum system. One 
may think of S as being the system of interest, like a particle enclosed in a container, 
and of C as a chain of measuring apparatuses £ that are brought into contact with the 
particle in a sequential manner. 

Our goal is to study the large time behaviour of repeated interaction quantum 
systems, and in particular, to describe the effect of the repeated interaction on the 
system S. One of our main results is the construction of the time-asymptotic state, 
which we call a repeated interaction asymptotic state (RIAS). 

States of S and £ are represented by vectors (or density matrices) in the Hilbert 
spaces Tls and Tie, respectively. We assume that dimTis < oo, while dimTCg < oo. The 
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observables of S and £ are bounded operators, they form the (von Neumann) algebras 
Tls C B{Hs) and StRg C BiJ-tg). Observables evolve according to the Heisenberg 
dynamics 1 1— > Tg{As) and 1 1— > rK^lg), where Tj and t| are groups of *automorphisms 
of SDTs and DJlg, respectively. 

We assume that there arc distinguished "reference" states, represented by the vec- 
tors ^1$ G lis and fif G 7^£, and for the purposes of the introduction, we shall take 
O5, to be equilibrium states with respect to Tj, r|, for inverse temperatures (5$, Ps, 
respectively. It is useful to pass to a description of the dynamics of vectors in Hs, He 
(Schrodinger dynamics). There are selfadjoint operators Lg, L£, called the standard 
Liouville operators, uniquely specified by 

r^(A) =e'*-^#y4e-"^#, and = 0, (1.1) 

where # stands here for either S or S. 

The Hilbert space of the entire system is given by H = Hs <8) He, where He, the 
Hilbert space of the chain, is the infinite tensor product ®rn>iHs- The non-interacting 
dynamics is defined on the algebra SJl^ <8'm>i by ®m>i T£- 

We consider interactions of the following kind. Fix an interaction time r > 0. 
During the interval [0,r), S interacts with the first element 6 in the chain C, while 
all other i?'s evolve freely. The interaction is specified by an operator V £ SJt^ Stt^. 
In the next time interval, [r, 2t), S interacts with the second element in the chain, 
through the same interaction operator V, and all other elements evolve freely, and so 
on. For t > we set 

t = m{t)T + s{t), (1.2) 

where m{t) is the integer measuring how many complete interactions of duration r have 
taken place at the instant t, and where < s{t) < r. We define the repeated interaction 
(Schrodinger) dynamics, for i > 0, ^ G H, by 

URi{t)tp = e-'<*)^m{t)+iQ-^rL^(t) . . . e-'-^^i^, (1.3) 

where 

is the generator of the total dynamics during the time interval [{m — 1)t, itlt). We have 
introduced L^, the operator on H that acts trivially on all elements of the chain except 
for the m-th one, and which, on the remaining part of H (which is just Hs ^ He), acts 
as 

L = Ls + Ls + V. (1.5) 

In (1.4), denotes the operator on H that acts nontrivially only on the A;-th element 
of the chain, on which it equals Lg. 

A state uj given by a density matrix on H is said to be normal. Our goal is to 
understand the time-asymptotics (i — > 00) of expectations 



u{Uni{trOUni{t)) =u{aij{0)), 



(1.6) 
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for normal states uj and for certain classes of "obscrvables" O. As mentioned above, we 
may regard S as the system of interest, so we certainly want to treat the case O G Tls. 
Another type of physical observable is of interest as well. Imagine we want to measure 
the variation, say, of the energy of 5 at a certain moment in time. This measuring 
process involves the system S, but also the element of the chain that is in contact 
with S at the given moment. We call such an observable an instantaneous observable. 
Various generalizations can be considered, see Section 2, but we limit our discussion in 
this introduction to the two kinds of observables just described. 

Asymptotic state. Let O be an instantaneous observable, determined by As 
and A^. This means that at time t = m{t)T + (1.2), O measures As ® A^ on 
the system S + £, where £ is the {m{t) + l)-th element in the chain C. We show in 
Theorem 2.3 that, under a natural assumption on the interaction, wc have 

u{a^m) - ^+{P<^thAs ® As)P) I 0, (1.7) 

as t — DO, where a;+ is a state on SJt^ which does not depend on oj (c.f. (1.13)), and 
where 

P = tns(^m>iPne, (1-8) 

with Pf2g denoting the orthogonal projection onto CQ,s- We identify the range of P 
with Hs- Relation (1.7) shows that the expectation of an instantaneous observable in 
any normal initial state approaches a r-periodic limit function (t i— > s(t) is r— periodic). 
The speed of convergence in (1.7) is exponential, ~ e"*'''/'^, where 7 > is a constant 
depending on the interaction. 

The restriction of the RIAS to the algebra of instantaneous observables character- 
ized by As G and As G SJtg is the r-periodic state 

As®As^ u+ {Pa^^ {As ® Ae)P) (1.9) 

on dJts®'^e-, see (1.7) (and (2.26) for the definition of the RIAS acting on more general 
observables) . 

The above-mentioned assumption on the interaction is an ergodicity assumption on 
the dynamics reduced to the system S. More precisely, we construct a (non— symmetric) 
operator K onH s.t. 

Pa\^i{As)Pns = {Pe'^^Pr^^^Pe'^^^^^AsPns, (1.10) 

with the property 

Pe'^^PVLs = O5. (1.11) 

We assume that 1 is a simple eigenvalue of Pe^'^^P, and that all the other eigenvalues 
lie strictly inside the complex unit disk. (We prove in Section 3 that this holds for 
concrete models). As a consequence of this assumption, we have 

(-pgirifp)m(t) ^ (1.12) 

as t — >■ 00, where JIJ is the unique vector in Hs satisfying {Pe^'^^ P)*fl^ = $7 J and 
{^g,^s) = 1- To arrive at result (1.7), where 

a;+(-) = {n*s, ■ ns) , (1.13) 
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we use (1.12) together with an argument involving a cycHcity property of Qg. 

Our approach is constructive in the sense that the asymptotic characteristics of the 
system, such as the state a;+, the speed of convergence 7, and the asymptotic dynamics 
in (1.7) can be calculated by rigorous perturbation theory (in V). 

Correlations & reconstruction of initial state. The (time dependent) asymp- 
totic expectations of observables do not depend on the initial state of the system, c.f. 
(1.7). However, asymptotic correlations do, and together with the asymptotic expec- 
tations they permit to reconstruct the intial state in the following way. 

Take a normal state u, an instantaneous observable O determined by As G SDts, 
As G 9Jl£, and an observable A G 971. We show in Theorem 2.5 that 



UJi 



{Aa'^,{0))-u;{A)u+{Pag>{As®A£)P)\^(}, (1.14) 



as i — > 00 (exponentially fast). According to (1.7) and (1.14), knowledge of the asymp- 
totic correlation function C+(t; A, O), and of the asymptotic expectation E+{t; O), de- 
termined respectively by 

lim \u{AaU{0)) - C+{t; A,0)\ = and lim \u;{ai-i{0)) - E+{t;0)\ = 0, 

t— >00 t—KX> 

allows for a reconstruction of the initial state uj according to 

Energy, entropy, average 2"*^ law of thermodynamics for RIAS. The 

formal quantity a^i{Lm{t)+i) , where is given by (1.4), has a well-defined variation 
in t. It is not hard to see by explicit calculation that this variation is zero in all time 
intervals [(m — 1)t, mr), and that it undergoes a jump 

j{m) = (Wl - Vrn) 

as time passes the moment mr. Here, denotes the operator V acting nontrivially 
only on Hs and the fc-th element Tig of the chain Hilbert space Tic- We interpret the 
variation of the above formal quantity as the (time dependent) observable of variation 
in total energy of the system. 

We show in Section 2.4.1 that for any normal initial state uj, the variation in energy 
during any time interval of length r takes the asymptotic expectation value a;+(j_|_), 
where 

j+ = PVP - Pal^j{V)P = -i r PaUlLs + Ls,V])Pds. (1.16) 

Jo 

(Here and in the rest of the paper we understand commutators to be defined in the 
form sense, but none of our arguments involve delicate domain questions with regards 
to commutators.) We define the (average) asymptotic energy production dE+ to be the 
change in energy during any interval of duration r, divided by r, in the limit of large 
times. This quantity is given by 

dE+ = -LO+iU) (1.17) 
r 
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and is independent of the initial state uj. We show in Section 2.4.2 that cj+(j+) > 0. 

Denote by loq the state on DJl determined by the vector Qs ®m>i & H, and let 
Ent(aj|u;o) denote the relative entropy of the normal state co w.r.t. loq. We think it 
is natural to define the entropy as a non-negative quantity, and our definition of it 
differs by a sign from the one given in [5]"^. We define the (average) asymptotic entropy 
production dS+ to be the change of (relative) entropy in any interval of duration r, 
divided by r, in the limit of large times. We prove in Section 2.4.2 that 

dS+ = ^u;+{U), (1.18) 
r 

where /?£ is the inverse temperature of the elements in the chain. The asymptotic 
entropy production does not depend on the initial state uj. We may combine (1.17) 
and (1.18) to arrive at an average 2"*^ law of thermodynamics for repeated interaction 
quantum systems, 

dE+ = T£dS+, (1.19) 

where Tg = 1/Pe is the temperature of the chain. Relation (1.19) is independent of the 
initial state of the system, and it holds for any repeated interaction system {V and r) . 

Repeated interaction quantum systems emerge in quantum optics, see e.g. [11, 4] 
and references therein. The effective evolution of the small system in certain specific 
regimes of parameters related to the Van Hove limit has been studied in [2]. In some 
parameter regimes repeated interaction models can be considered as coarse grained 
versions of a system S in contact with quantum noises [3] . 



2 Model Results 

2.1 Repeated interaction models 

The models we consider consist of a system S which is coupled to a chain C = £+£ + ■■ ■ 
of identical elements £. We describe S and £ as 14/'*-dynamical systems {dJls,T^) and 
(9Jtf,r£-), where Tls, ^fftg are von Neumann algebras "of observables" acting on the 
Hilbert spaces Hs, He, respectively, and where and r| are (cr-weakly continuous) 
groups of * automorphisms describing the Heisenberg dynamics. In this paper, we 
consider the situation dimTCg < oo and dimTCg < oo. 

We assume that there are distinguished vectors Qs £ ^5 and Of G Tig, determining 
states on dJls and DJlg which are invariant w.r.t. Tj and t|, respectively, and we assume 
that Qs and Qg are cyclic and separating for Tls and Tig, respectively. One may 
typically think of these distinguished vectors as being KMS vectors. 

The Hilbert space of the chain C is defined to be the infinite tensor product 

nc = ®m>ing (2.1) 

w.r.t. the reference vector 

Vtc = ^g®0.g-- - . (2.2) 

^For a finite system we have Ent(a;|a;o) = Tr(p(log p — log po)), where p and po are density matrices 
determining the states u! and wo, respectively, and where po > 0. 
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In other words, Tic is obtained by taking the completion of the vector space of finite 
linear combinations of the form 0m>i'4'm, where ■i/'m £ T^£, i^m = except for finitely 
many indices, in the norm induced by the inner product 

{®mi^m, ®mXm) = Hm (V'm, Xm)^^ • (2-3) 

We introduce the von Neumann algebra 

SWc = ®m>me (2.4) 

acting on (8'm>i'^£) which is obtained by taking the weak closure of finite linear com- 
binations of operators (8)m>i^m5 where Am G DJlg and Am = except for finitely 

many indices. 

The operator algebra containing the observables of the total system is the von 
Neumann algebra 

m = ms^mc (2.5) 

which acts on the Hilbert space 

n = ns^nc. (2.6) 

The repeated interaction dynamics of observables in dJl is characterized by an inter- 
action time < r < 00 and a selfadjoint interaction operator 

VeTls^me. (2.7) 

For times t € [r(m — l),rm), where m > I, S interacts with the m-th clement of the 
chain, while all other elements of the chain evolve freely (each one according to the 
dynamics rg). The interaction of S with every element in the chain is the same (given 
by V). 

Let Ls and be the standard Liouville operators ("positive temperature Hamil- 
tonians", c.f. references of [6, 9]), uniquely characterized by the following properties: 
(where # = S,£) are selfadjoint operators on H# which implement the dynamics 

4(^4) = e'*^#y4e-'*^#, e (2.8) 

and 

L#n# = 0. (2.9) 

We define the selfadjoint operator 

L = Ls + L£ + V, (2.10) 

omitting trivial factors Is or Is (by Ls in (2.10) we really mean Ls (8) etc). L 
generates the automorphism group e'*^ ■ c^'*^ of dJls the interacting dynamics 
between S and an clement £ of the chain C. The explicit form of the operator V is 
dictated by the underlying physics, we give some examples in Section 3. 
For m > 1 let us denote by 

Lm = Lm+ ^ Ls^k (2-11) 
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the generator of the total dynamics during the interval [(m — l)r, mr). We have intro- 
duced Lm, the operator on H that acts trivially on all elements of the chain except for 
the m-th one. On the remaining part of the space (which is isomorphic to Hs <8) Hs), 
Ljn acts as L, (2.10). We have also set Lg^fc to be the operator on H that acts non- 
trivially only on the A;-th element of the chain, on which it equals Lg . Of course, the 
infinite sum in (2.11) must be interpreted in the strong sense on Ti. 
Decompose t G as 

t = m{t)T + s{t), (2.12) 

where m{t) is the integer measuring the number of complete interactions of duration 
r the system S has undergone at time t, and where < s(t) < t. The repeated 
interaction dynamics of an operator A on H is defined by 

aii{A) = Uni{t)*AUni{t) (2.13) 

where 

Ujii{t) = e-'<*')^m{t)+iQ-^rL^(t) . . . Q-irLi ^2.14) 

defines the Schrodinger dynamics on H. According to this dynamics S interacts in 
succession, for a fixed duration r and a fixed interaction V, with the first ■m{t) elements 
of the chain, and for the remaining duration ,s{t) with the {m{t) + l)-th element of the 
chain. Being the propagator of a "timc-dcpcndcnt Hamiltonian" (which is piecewise 
constant), C/Ri(t) does not have the group property in t. 

Our goal is to examine the large time behaviour of expectation values of certain 
observables in normal states lo on dJl (states given by a density matrix on H). The 
system S feels an effective dynamics induced by the interaction with the chain C. Under 
a suitable ergodicity assumption on this effective dynamics the small system is driven 
to an asymptotic state, as time increases. We will express the effective dynamics and 
the ergodic assumption using the modular data of the pair {Tls 9Jtc, 0,s <Si^c)- 

Let J and A denote the modular conjugation and the modular operator associated 
to (arts (g) SJtf, (g) Qf), [5]. We assume that 

(A) AV2yA-V2 e $H£ 

and we introduce the operator 

K = L- JA^/VA-^/^ J, (2.15) 

called a C-Liouville operator, [6, 9]. It generates a strongly continuous group of 
bounded operators, denoted e'*^, satisfying ||e'*^|| < el*l H^^''^^^ ^''^H. The main fea- 
ture of the operator K is that c^^^ implements the same dynamics as e'*^ on Stts (g QJlg 
(since the difference K — L belongs to the commutant (8) SJig), and that 

K^S ®^£ = Q. (2.16) 

Relation (2.16) follows from assumption (A), definition (2.15) and the properties A~^/^ J = 
JA^/^ and JA^I'^A^s ®^s = A*^s ® ^S, for any A&ms®Ms- 
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Let 

P = tns(^\^c){^c\ (2.17) 

be the orthogonal projection onto TCs (g) CQc — 'Hsi where Qc is given in (2.2). If B is 
an operator acting on TC then we identify PBP as an operator acting on 7is- We have 

Proposition 2.1 There is a constant C < oo s.t. \\{Pe^^^ P)'^\\b{'Hs) — '^^^ 
t G M, m > 0. In particular, spec(Pe^*^P) C G C | < 1} and all eigenvalues 
lying on the unit circle are semisimple. 

We give a proof of Proposition 2.1 in Section 4.3. Relation (2.16) implies that for all 
i G M, Pe^^^P^s = ^S- Our assumption (E) on the effectiveness of the coupling is an 
ergodicity assumption on the discrete dynamics generated by 

M = M{t) = Pq^'^^P (2.18) 

(E) The spectrum of M on the complex unit circle consists of the single eigenvalue 
{1}. This eigenvalue is simple (with corresponding eigenvector ^s)- 

Assumption (E) guarantees that the adjoint operator M* has a unique invariant vector, 
called fig (normalized as (0^,05) = 1), and that 

lim M"* = 7r:=|J^5)(J^5|, (2.19) 

m— >c» 

in the operator sense, where vr is the rank one projection which projects onto CVts 
along (COJ)^. In fact, we have the following easy estimate (valid for any matrix M 
with spectrum inside the unit disk and satisfying (E)) 

Proposition 2.2 For any e > there exists a constant C\ s.t. HM™— 7r|| < Cee~"^^"'~'^\ 
for all m>0, where 7 := min_j,gspgc(M)\{i} I log kl I > 0. 

The parameter 7 measures the speed of convergence. 

Remark. If all eigenvalues of M are semisimple then in Proposition 2.2 we have 
||M"^ — 7r|| < Ce~™^ for some constant C and all m > 0. 

As a last preparation towards an understanding of our results we discuss the kinds 
of observables we consider. One interesting such class is QJlg C 9Jl which consists of 
observables of the system S only. There are other observables of interest. We may 
think of the system S as being fixed in space and of the chain as passing by S so that 
at the moment t, the {m{t) + l)-th element £ is located near S, c.f. (2.12). A detector 
placed in the vicinity of S can measure at this moment in time observables of S and 
those of the {m{t) + l)-th element in the chain, i.e., an "instantaneous observable" of 
the form ^4^ (g) '&m{t)+i{^o), where As G SDts, Bq G SDtf, and -dm : 3Jtf — SDTc is defined 

by 

^miAs) = l£---l£<S)A£0ls--- (2.20) 

where the Ag on the right side of (2.20) acts on the m-th factor in the chain. An 
example of such an observable is the energy flux (variation) of the system S. More 
generally we may be interested in the expectation value of operators of the form 



[As; Ai; Bj] = Ag <=i Ai 'dm{t)+j+i{Bj), 



(2.21) 
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where As G Tig, Ai, . . . ,Ap £ Tig, B_i, . . . , Bq, . . . , Br G Tig and where t = m{t)T + 
sit) as in (2.12) and 'dk is given in (2.20). The parameters p > 1, r > are not 
displayed in the l.h.s. in (2.21). (We always assume that p < m{t) — £ + 1.) As and 
the Ai represent observables we measure on the small system and on the element with 
index i of the chain, the Bq is the "instantaneous" observable, measured in the element 
m{t) + 1 of the chain (the one in contact with S at time t), while the Bj with negative 
and positive index are the quantities measured in the elements preceding and following 
the {m{t) + l)-th. 

2.2 Asymptotic state 

Throughout the paper we assume that Conditions (A) and (E) of the previous section 
are satisfied. 

We consider the large time limit of expectations 

E{t)=uj{a'^^[As;Ai-Bj]) (2.22) 

for observables [As; Af, Bj] as in (2.21) and for normal states u on 93t. Define the state 
a;+ on DJts by 

u;+{As) = {n*s,Asns), (2.23) 

where fij is defined in (2.19). 

Theorem 2.3 Let lo be fixed and take A^ = Ig, i = 1, . . . ,p. For any e > there is a 
constant Cg s.t. for all t >0 

\E{t) - E+{t)\ < Cee-*(^-^)/^, (2.24) 

where j > is given in Proposition 2.2, and where £'+ is the t -periodic function 

E+it) = C0+ (PaS"+^(*) (^5 ® ® • • • ® Bo) P) nU (Bj)^^ • (2.25) 
Here, {B,)^^ = {^s, Bj^g) . 

We define the RIAS to be the r-periodic state on Tts <^j=_£ ^£ given by 

As Bj ^ u+ (Pa5[+^(*) {As®B_t®---® Bo) p) UU {Bj)n, ■ (2-26) 

Using (2.24) and the unicity of the limit, one can see that actually the state a;+ does 
not depend on the choice of the reference state O5. 

Remarks. 1) If B^i, . . . , -B_^'_i = for some —i' — 1 < —1, then one shows that 

iu+ ® ® • • • ® Bo) P) = a;+ [Pa^l^'^''^ (^Ag <^ ■ ■ ■ <^ Bo)p) , 

£T-\-s(t) 

and in case Bj = for all j = —£, . . . , formula (2.25) is understood with ccj^j 
replaced by a^l^ . 
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2) Ce in Theorem 2.3 is uniform in r for r > varying in compact sets, and it is 
uniform in {As G Ms, {BjYj^^ C Ms \ \\As\\ irj=i\\Bj\\ < const.}. 

3) The convergence is determined by that of Proposition 2.2. If the ergodic as- 
sumption (E) is not satisfied then the hmit lim„_>oo M" still exists, in a weaker sense. 
Namely, if there arc eigenvalues different from 1 on the circle, then the limit exists 
in the ergodic mean sense, jf^n=o^^ ~ ^ + ^(ir)- Further, if 1 is a degenerate 
eigenvalue of M then the limit exists but the projection ir is not one dimensional. This 
reflects in Theorem 2.3 in the following way. If 1 is non degenerate, but there are other 
eigenvalues on the circle, then Theorem 2.3 holds with (2.24) replaced by 



t ^ — \ T 



m=0 



< y- (2.27) 



If on the other hand 1 is degenerate but there is no other eigenvalue on the circle, 
then one can still prove that the expectation value E{t) has an aymptotic behaviour 
Eoo{t,Lo), which is r— periodic, but which will a priori depend on the initial state to 
(c.f. (4.25) in the proof of Theorem 2.3). Of course if both 1 is degenerate and there 
are other eigenvalues on the circle, then one gets convergence to E^{t,uj) but in the 
ergodic mean. 

Our next result incorporates the measurement of obscrvables Ai, . . . Ap G 971^ for 
a chain consisting of dispersive systems £. We measure dispersivity by the property of 
return to equilibrium. £ is said to have the latter property iff for any normal state lo£ 
on SDtg we have the relation 

lim LosiT'eiAs)) = (Qf , AeQe) , (2.28) 

t— >oo 

for any Ag G 971^. Examples of such £ include reservoirs of ideal quantum gases. It 
is worthwile to mention that £ has the property of return to equilibrium if and only if 
e'*^^ converges in the weak sense to the orthogal projection onto CQ^, as i — > oo. 

Theorem 2.4 Suppose £ has the property of return to equilibrium. Then 

lim I £;(*)- £;+(*)! =0, (2.29) 

t— >oo 

where E^(t) is the t -periodic function 

E+{t) = a;+(Pag"+^(*) (^5 <8) ® • • • ® ^0)^)^=1 (^i)Q, UU (^^K ■ (2-30) 

Remark. The speed of convergence in (2.29) is determined by that of return to equi- 
librium, (2.28), and by 7, Proposition 2.2. The limit (2.29) is uniform in r, for r varying 
in compact sets, and it is uniform in balls of observables \\As\\ nf=i IIA|| 11^=-^ W^jW ^ 
const. 
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2.3 Correlations &; reconstruction of initial state 

As Theorems 2.3 and 2.4 show, the hmiting expectation values E^(t) are independent 
of the initial state (the state is, c.f. (2.23)). However, limiting correlations are not, 
and their knowledge allows to reconstruct the initial state. 

Fix a normal initial state w of 93t and let A G SDT, As G ^XHs, Bq G VJle. We define 
the correlation between A and the instantaneous observable As '!?m(t)+i(-^o) by 

C{t;A,As,Bo)=u;{Aal,,{As^^n.it)+i{Bo))). (2.31) 

Theorem 2.5 For any e > there is a constant Cg s.t. for all t > 

\C{t;A,As,Bo)-C+it;A,As,Bo)\ < Qe-*M/r^ (2.32) 

where 7 is given in Proposition 2.2, and where C+ is the t -periodic limiting correlation 
function 

C+{t- A, As, Bo) = io{A) u+ (Paf^ {As 5o)p) , (2.33) 

with a;+ defined in (2.19). 

Remark. Relation (2.33) allows us to reconstruct the initial state w, knowing the 
asymptotic state a;+ and the asymptotic correlation function C+. 

2.4 Energy, entropy and their relation 

It may not be meaningful to speak about the total energy of the system, because it 
may have to be considered as being infinite, e.g. if the elements £ of the chain are 
infinitely extended quantum systems with non-vanishing energy density. However, we 
can define the time variation of the total energy of the system and link it to its entropy 
variation, giving us an average 2"*^ law of thermodynamics for RIAS. 

2.4.1 Energy 

Recall that L^+i, m > 0, is the generator of the total dynamics in the time interval 
t = rnr + s G [mr, (m + 1)t), during which the (m+ l)-th element of the chain interacts 
with 5, c.f. (2.11). Given any integer m > and any < s < r it is easy to formally 
verify the relation 

aS"+*-(Wi) - «Rf (^m+i) = 0. (2.34) 

This suggests that the formal quantity aj^j(Lf„_|_i) is constant for t in any interval 
[mr, (m + 1)t). Another short calculation yields that this quantity undergoes a jump 
j{k) as time passes the moment kr, k>l: for {k — l)r <t\<kT<t2<{k + l)r we 
have 

j{k) := a'i,{L,+,) - a^i(Lfc) = a'^,{Vk+i - V^), (2.35) 

where we set 

Vk = [lms^^k]{V) (2.36) 

(see (2.20)). We interpret j{k) as the change in total energy as time passes the moment 
kr. 
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Theorem 2.3 tells us that for any normal state a; on 9Jl and for any e > 0, there is 
a constant s.t. 

|u;(j(fe))-a;+(j+)|<C,e-^(^-), (2.37) 

where 

U = PVP - Pa-^,{V)P = -i r Pal,,{[Ls + L£,V])Pds. (2.38) 

Jo 

Relation (2.37) and the fact that the energy is piecewise constant shows that uj^{j-^-) 
is the change of energy in any interval of length r, in the large time limit. We thus call 

dE+ = -io+iU) (2.39) 
r 

the asymptotic energy production. The asymptotic energy production does not depend 
on the initial state of the system. 

Remark. It is not hard to see that the expectation of the energy jump is constant 
in the state L0-^- loc, where uc is the vector state on dJlc determined by ^Ic, (2.2): 

u;+ ® ucUik)) = uj+ij+), VA; > 1. (2.40) 

We introduce the variation of the total energy, AE{t), between the instants t = 
m{t)T + s{t) and t = 0. It is the sum of the energy jumps, 

m{t) 

AE{t) = J2j{k), iovt>T, (2.41) 
k=l 

and AE{t) = if < t < r. Estimate (2.37) shows that for any normal state a; on 93t 
there is a constant C s.t. 

t 

for all t > 0. The energy grows asymptotically linearly in time. 
2.4.2 Entropy, average 2"^ law of thermodynamics 

Let to and cuq be two normal states on dJt. The relative entropy of uj with respect to (Jq 
is denoted by Ent(u;|a;o); where our definition of relative entropy differs from that one 
given in [5] by a sign, so that in our case, Ent(a;|a;o) > 0. 

For a thermodynamic interpretation of the entropy and its relation to the energy 
variation, we assume in this section that 0,^ is a (/35,r^)-KMS state on SlJt^, and that 

is a (/3£,r|)-KMS state on Tie, where Ps^Ps are inverse temperatures. Let loq be 
the state on 9Jl determined by the vector (g) Qc (c-f- before (2.1), and (2.2)). 

We are interested in the change of relative entropy of the repeated interaction 
system as time evolves. 

Proposition 2.6 Let lu be any normal state on 9Jt. Then Ent(a;oQ;j^j|a;o) is a contin- 
uous, piecewise differentiable function oft>0. Moreover, we have 

Ent(a; o a^ilwo) — Ent(a;|a;o) = w {psAE{t) - a^i(X(i)) + X(0)) , (2.43) 



< 



C 



(2.42) 
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where AE{t) is the variation of the total energy between the moments t = and t = 
m{t)T + s{t), (2.41), and where 

X{t) = P£V^(^t)+i + iPs - Ps)Ls, (2.44) 

with Vfc given by (2.36). 

The proof of (2.43) is based on the entropy production formula [7]. We give it in Section 
4.4. It is not hard to verify that for t G (mr, {m + l)r) we have 

_d 

dr 

and that left and right derivatives of Ent(a; o Q;^j|a;o) exist as t — mr, but they do not 
coincide. 

If Ent(a;|ti;o) < 00 then all terms in (2.43) are bounded uniformly in i, except 
possibly Ent(w o a^jjwo) and io{f3sAE{t)). Hence (2.39) and (2.42) show that for any 
normal state a; on SDt there is a constant C s.t. 



-Ent(a; o aRi|a;o) = -u{a^i(\[(3sLs + /^fl'f ,m(t)+i, Hn(t)+i])) , (2-45) 



Ent(a;oaRT|a;o) 

^+{3^ 



< J, (2.46) 



t T 

for all t > 0. The entropy grows linearly in time, for large times. 
The relative entropy is non-negative, so (2.46) shows that 

^+(i+) > 0. (2.47) 

We show in Section 3 that a;+(j+) is strictly positive for concrete systems. It follows 
from (2.46) also that 

sup |Ent(a; o or^jIwo)! < 00 '^+(j+) = 0- (2.48) 

Since a;+(j+) is independent of lo it follows that for a given interaction (y, r) the relative 
entropy either diverges for all initial states a;, as i — >■ 00, or it stays bounded for all 
initial states u. In particular, i/ti;+(j+) > then there does not exist any normal state 
a; on 9Jl which is invariant under a\^i (i.e., such that lo o = uj, for all t > 0). 

Proposition 2.7 We have 

lim [Ent{u o a*^/\LOo) - Ent(a; o a^ilwo)] = /3£ w+(j+)- (2.49) 

The change of entropy during an interval of duration r, for t ^ 00, is thus given 
by /?£ u+{j+) > 0. We call 

dS+ = ^u+iU) (2.50) 
r 

the (average) asymptotic entropy production. The quantity dS-^- represents the increase 
in entropy per unit time, in the limit of large times. It does not depend on the initial 

state of the system. 

Remark. One sees easily that the expectation of dS+ is constant in the state 
a;+ uoc (see also (2.40)). 

Relations (2.50) and (2.39) lead us to the average 2^'^ law of thermodynamics, 

dE+ = T£dS+, T£ = l/Pe. (2.51) 

This law does not depend on the initial state of the system. 



BJM February 7, 2008 



14 



3 Examples 

3.1 Spin-Fermion system with quadratic interaction 

As our main example, we consider the case where the small system <S is a 2— level system 
and the elements of the chain consist of free Fermi reservoirs at positive temperature 
(3~^. Let us first describe precisely the model (see also [6] and references therein). 
The von Neumann algebra of observables for the small system is 

Ms = M2(C) ^l = {A<^l\Ae M2(C)} (3.1) 

acting on the Hilbert space 

7^5 = (8) C^. (3.2) 
Let ax, (Jy, Oz be the usual Pauli matrices, i.e. ~ ^\^\ 0)'*^^"^ ^i o)' 
= ^ ^ ^1 ) '^'^^ dynamics of the small system is then given by 

T% {A®t) = ^^"^ Ae-'*"^ ® 1. (3.3) 

For convenience we chose the reference state 0^5 to be the tracial state, i.e. usiAiSi 1) = 
^Tr(^). Note that it is a (r5,0)— KMS state. Its representative vector is 

1 1 
ns = -j=ipi (8) V'l + -^i^2 <8) ^2 (3.4) 

where {'ipi,ip2) is the canonical basis of C^. For shortness, we will denote by ipij := 
V'i ® ipj the corresponding basis of Hs- The standard Liouvillean then writes 

Ls = (Tz®l-l0az, (3.5) 

and its spectrum is spec(L5) = {—2,0,2} where has multiplicity 2, and —2, 2 are 
non degenerate. Finally, the modular conjugation and modular operator associated to 
{Ms, ^s) are 

Js{(j)®'^) = il^®(j), A5 = 1(8)1, (3.6) 

and where " denotes the usual complex conjugation on C^. 

We then describe an element of the chain, i.e. a free Fermi gas at inverse temper- 
ature (3. Let [) be the Hilbert space of one single fermion and h its energy operator. 
The operators a(/) and a*(/) denote the corresponding annihilation and creation op- 
erators acting on the fermionic Fock space r_([)) and they satisfy the canonical anti- 
commutation relations (CAR). As a consequence of the CAR, the operators a(/) and 
a*{f) are bounded and satisfy ||a'^(/)|| = ||/|| where a* stands either for a or for 
a*. The algebra of observables of a free Fermi gas is the C*-algebra of operators 21 
generated by {a*(/)|/ G f)}. The dynamics is then given by 

T/(a#(/)) = a#(e^*V)- (3-7) 
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It is well known (see e.g. [5, 10]) that for any /3 > 0, there is a unique (rf , KMS 
state LUg on 21 which is determined by the two point function a;^(a*(/)a(/)) = (/, (1 + 
g/3ft)-iy^ Finally, let be the Fock vacuum and N the number operator. 

We now fix a complex conjugation (anti-unitary involution) / — ^ / on [) which 
commutes with the energy operator h. It naturally extends to a complex conjugation 
on the Fock space r_(f)) and we denote it by the same symbol, i.e. $ ^ 

The GNS representation of the algebra 21 associated to the KMS-state LOf) is the 
triple {Hf,t^/3,^f) [1] where 

= r_(f)) ®r_(i)), OF = f^f0f^f, (3.8) 

and 



(3.9) 



The von Neumann algebra of observables for an element of the chain will then be 
the enveloping von Neumann algebra 

m = 7r/3(2l)", (3.10) 

acting on the Hilbert space 

He = Hf. (3.11) 

The dynamics on 7r^(2t) is given by 

rUMA)) = MrkA)) (3.12) 

and extends to SJtf in a unique way. The representative vector of the equilibrium state 
is 

ns = Qf, (3.13) 
and the standard Liouvillean then writes 

= dr(/i) (g) 1- ]i®dr(/i). (3.14) 

Finally the modular conjugation and the modular operator associated to (Tlsjils) are 
j£($®*) = (-l)^(^-i)/2^<^(_l)^(^-i)/2$^ A£ = e-^-^^. (3.15) 

We finally specify the interaction between the small system and the elements of the 
chain, i.e. the operator V. Let g e t) he a form factor, we set 

V :=a^(g) lc2 ® apig)af^{g) e (3.16) 

This is the simplest non trivial interaction for which the number of particles is con- 
served. 

We moreover assume that 
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(SFl) e'^'^/^g e f). 

This ensures that Assumption (A) is satisfied. Indeed, using (3.6)-(3.15)-(3.16), we get 

^1/2^^-1/2 ^3_;L7) 



1 \ . , ..AT- / e-«V2 



g] 01+ (-1)^ ® a , 



X a a 1+ (-1)^ ® a* ' 



]1C2 



The Liouville operator which generates the interacting dynamics is then the selfadjoint 
operator 

Lx := Ls + L£ + XV, (3.18) 
while the C-Liouville operator is 

Kx := Ls + L£ + X{V - JA^/VA-V2 J) = Kg + XW, (3.19) 

and where A G M is a couphng constant. 

For simphcity reasons we will moreover assume that 

(SF2) [) = L2(R+,g) where g is some auxiliary Hilbert space and the operator h is the 
multiplication operator by r G M+. 

Finally, let gpir) := (1 + e-^'')-V2^(r). 

Note also that the system S has the property of return to equilibrium. 
Our first result is the 

Theorem 3.1 Suppose Assumption (SF1)-(SF2) are satisfied. Then for any t ^ f + 
ttN, there exists Aq > such that for all < |A| < Aq, the operator Mx := Pe^'^^^P sat- 
isfies the ergodic assumption (E). In particular the spin-fermion system with quadratic 
interaction satisfies Theorem 2.3, with 7 = ^ ("i^+"2) ^2 _|_ (9(A'^), Theorem 2.4, and 
moreover the asymptotic state lo+^\ is given by 

uj+x{As) = 7 (aiV'ii + a2'i/'22,^5(V'ii +'^22)) (3.20) 

ai + a2 



+A 

where, for j = 1,2, 



Q-ph/2 



5/3 



t) 2{ai + 02 



■(V'12 + ^^21, ^5(V'll + V'12)) + O(A^), 



aj :=JJ dridr2e-'^^^||5^(ri)||2||5^(r2)||2sinc2 (^l(?_^±I^j > q, (3.21) 
sinc(a;) := and all the integrals run overR.'^. 
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The unperturbed operator Mq has eigenvalues 1 (with multiplicity 2), e^'"^ and e~^'^ 
(see (3.5)). The assumption on the interaction time r ensures that these eigenvalues do 
not coincide and makes the computation in perturbation theory as simple as possible. 
However, it can probably be omitted. 

One can also see that the asymptotic expectation E^(t) (see (2.25)) has a non 
trivial periodicity (i.e. it is not constant) at the order A. 

We now turn to the question of entropy production for this simple model. We will 
prove that it is strictly positive, at least for small coupling constant. More precisely, 
we have 

Theorem 3.2 Suppose Assumptions (SF1)-(SF2) are satisfied. Then for any r ^ 
^ + ttN, there exists Ai > such that, for all < |A| < Ai, the system has strictly 
positive asymptotic entropy production. 

3.2 Spin-Fermion system with linear interaction 

As a second example, we consider the same spin fermion system at inverse temperature 
p. The only change concerns the interaction term. Let ^ G f) be a form factor, we set 

V:=a^® lc2 ® {ap{g) + a}{g)). (3.22) 

Once again, we assume that Assumption (SFl) holds which ensures that G 
OJl. Indeed 

^1/2^^-1/2 

We then have the same kind of result as for the case of quadratic interaction, namely 

Theorem 3.3 Suppose Assumption (SF1)-(SF2) are satisfied. Then for any r ^ 
I + 7rN, there exists Aq > such that for allO < |A| < Aq, the operator Mx := Pe^'^^^P 
satisfies the ergodic assumption (E). In particular the spin-fermion system with linear 
interaction satisfies Theorem 2.3, with 7 = ("1^+02) ^2 _|_ 0{\^), Theorem 2.4 and the 
asymptotic state co+^x is given by 

^+,xiAs) = ] (aiV-ii + ^21^22, Asii^n + ^22)) + ©(A^), (3.23) 




where 

ai := J dr\\g/3{r)\\l (^e-^^'sinc^ (li^^^^j + sine 
"2 := J dr||5^(r)||^ ^^e-'^^'sinc^ ^ r(r + 2) ^ _^ 

and all the integrals run over M"*" . Moreover the system has strictly positive asymptotic 
entropy production. 
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3.3 Spin-Spin model 

As our last example, we consider a model in which the small system as well as the 
elements of the chain consist of a 2— level system. Such kind of systems (or more 
generally a d— level system interacting with a chain of n— level systems) have been 
considered previously in [2] . 

The von Neumann algebra of observables for the small system and for the elements 
of the chain is 

ms = ms = M2(C) ^1 = {A(^l\Ae M2(C)} (3.24) 
acting on the Hilbert space 

ns = ne = c'^® (3.25) 

Let Es and E£ be non negative real numbers. They will play the role of the energy of 
the "excited" state of the small system and of the elements of the chain respectively. 
The dynamics of the small system is then given by 

T^{A ® 1) = e'^^^Ae-'^''^ ® H, (3.26) 

and the one of an element of the chain by 

4{A (g) 1) = e'^^^ Ae-'*''^ 1, (3.27) 

where/.5=(Q E^)^^dhe=(^l 

Once again, for convenience, we chose the reference state lo^ to be the tracial state, 
i.e. uJs{A (g) 1) = iTr(yl) (the results of course do not depend on this choice). Its 
representative vector is (Section 3.1) 

= ^V'll + ^^22. (3.28) 

The standard Liouvillean writes 

Ls = hs<^l-l<^hs, (3.29) 

and the modular conjugation and modular operator associated to {Tls, ^s) are 

Js{(t>®il^) = ^®4>, A5 = 1(8)11. (3.30) 

In order to avoid confusions between the small system and an element of the chain 
we will denote by (t)ij = 4>i® 4>j instead of V'ij the basis of He . The reference state us 
will be the (t£,/3)— KMS state. Its representative vector writes 

= ^ (0n + e-^^^/2</,22). (3.31) 
V 1 + e-l^^£ 

The standard Liouville operator is 



(3.32) 
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and the modular conjugation and modular operator associated to {WlsjQs) are 

Je{(f>^tp) = tp^^, Af = e-^-^^ . (3.33) 

Note that here the system £ does not have the property of return to equilibrium. 

We now describe the interaction. Let us denote by a and a* the annihilation and 
creation operators associated to the vectors 4>i (ground state) and 02 (excited state), 

i.e. a^i = 0,a^2 = (t>i,a*(l)i = (l)2,a*(f>2 = 0. Finally let J = ^ ^ d ) ^ -^2(C). The 
interaction is then given by 

V = I^l^a*®l + r^l^a®l. (3.34) 

The Liouville operator which generates the interacting dynamics is then the selfadjoint 
operator 

La := Ls + L£ + XV, (3.35) 
while the C-Liouville operator is 

Kx := Ls + L£ + X{V - JAV2yA-V2 J) = Ko + XW, (3.36) 

where A is a coupling constant. 

We finally consider the following assumptions. 

(551) b 7^ and t{E£ - Es) ^ 27rZ. 

(552) c 7^ and T{Ee + Es) ^ I'k'L. 

If either (SSI) or (SS2) is satisfied, then the ergodic assumption (E) holds. 

Theorem 3.4 Suppose that rEg ^ ttZ and that either Assumption (SSI) or (SS2) 

is satisfied. Then, there exists Aq > such that for all < |A| < Aq, the operator 
M\ := Pe^^^^P satisfies the ergodic assumption (E). In particular the spin-spin system 
satisfy Theorem 2.3, with 7 = 7oA2 + 0(A'^), and the asymptotic state w+,a given by 

Lo+,xiAs) = ] (aiV'ii + 02^^22, Asii^u + ^^22)) + O(A^), (3.37) 

ai + 02 



where 

0.2 



V2(ai + a2) T2(a]^ + a2) 

7o := mm ' 



1 + e-/3^£ ' 2(1 + e- 



+^sinc2(l|£)(|a|2 + |d|2 - ad - ad) 

// moreover both Assumptions (SSI) and (SS2) are satisfied, then the system has 
strictly positive asymptotic entropy production. 
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Once again, we assume that tE^ ^ vrZ just in order to make the eigenvalues of 
Mq not coincide and this can probably be weakened. On the other hand, Assumptions 
(SSI) and (SS2) are much deeper. Their signification is that there is an effective 
coupling between the ground state and the excited state of the small system (b or c non 
zero) as well as a non resonant phenomenon between the energies of the small system 
and the elements of the chain {t{E£ — Eg) or t{E£ + Es) not in 27rZ). Asking that 
either (SSI) or (SS2) is satisfied is actually equivalent to the condition ai + a2 / 0. 

4 Proofs 

4.1 Proofs of Theorems 2.3, 2.4 

We give the full proof of Theorem 2.4, the proof of Theorem 2.3 is a special case of the 
former. 

It is enough to show (2.29), (2.30) for vector states = ■ ifj), ip e H, = 
1. Further, since every ip £ Ti. is approximated in the norm of Ti by finite linear 
combinations of vectors of the form ipg <X'm>i i^rm where ^5 £ TCg, = ii m > N , 
for some N < 00, it suffices to prove (2.29), (2.30) for vector states determined by 
vectors of the form 

V'5 «)m=l V-m ®m>N € H, (4.1) 

where = ||^m|| = 1, 1 < m < iV, for arbitrary N < oo. Finally, since the 

vectors Us, are cyclic for the commutants 971^, SDt^, any vector of the form (4.1) is 
approximated by a 

il; = B' fls^ nc, (4.2) 

for some 

B' = B's0Z=lBL^^rn>NT^£ ^ Tt' , (4.3) 

with B'g e 9Jl^, B'^ 6 m'^. It is therefore sufiicient to show (2.29), (2.30) for vectors 
of the form (4.2), (4.3). 

Let As G 3Jt5, Ai,... ,Ap G dJl^ and B^i, . . . , Bq, . . . ,Br & dJls be fixed observables 
{i,r > 0). We examine the expectation value 

m ■■= {A «RI (^5 ®Li Ai l}m(t)+J+l{Bj)) V> , (4.4) 

where V' is given in by (4.2), m{t) is determined by (2.12), and where the Ai act on 
the first p factors of TCq. It is clear that the trivial factors l^- are omitted in (4.4). The 
choice of indices in (4.4) is such that at time t the observable Bq is measured in the 
element 6 of the chain which is interacting with S (i.e., the (m(t) + l)-th element of 
the chain). The following decomposition serves to isolate the dynamics of the elements 

5 which do not interact at a given time, see also (2.11). 

g-isWig-irL^ . • • e"'""^! = C/-e"^*-^"'+ie~"'''^'" • • • e~'''^^U^, (4.5) 
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where 



Um = exp 



= exp 



m+l 

-i ^ (j - l)TLf J - i(mr + s) ^ Lsj 

3=2 j>m+2 



(4.6) 
(4.7) 



We obtain from (2.13), (2.14) and (4.5) 



E{t) = {i), {U^ y e 



m ) 



where we write m,s for m{t),s{t). The operator {U^)* ■ ■ - U^ in the r.h.s. of (4.8) 
belongs to 9Jt. Hence we can commute the B' in (4.2) with this operator to obtain 



E{t) = (O, {B')*B'{U+)*e'^^^ ■ ■ ■ 



(m— i)r+s 



y g— isL^+i g— irLn 



where we have set 



e---^^n)fl{B,)^^ 



{-j-l)T+S 



{Bj))^l9rn+liBo) 



(4.9) 



(4.10) 

and we write (O)^ = (x, Ox) for an operator O and a vector x- The operator 
disappears because it leaves Q invariant, c.f. (2.9). We are able to factorize the averages 
(Bj)^^, for j > 1, because B' and all propagators in (4.8) act trivially on factors of TCc 
with index > m + 2 (note also that N < m since we have in mind the limit m oo). 

Since {B')*B' acts trivially on factors in He with index > iV + 1 we may replace 
{U^r = (C^i))* in (4.9) by 



Un = exp 



N 



(4.11) 



which is a unitary not depending on t. Without changing the value of (4.9) we can 
replace in that equation successively Li by Ki, then L2 by K2, up to replacing Lp by 



Kp. Then 



E{t) = {n, {B'YB'Un e'^^^ ■ ■ ■ e'^Kpe'^^v+^ . . . ^^rLm^isLm+i 



,-1 
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where we have used that 17 is in the kernel of the Kj, and where wc have commuted 
the product of the freely evolved Aj^s to the right through the propagators which act 
on different factors. 

Because the system S has the property of return to equilibrium the propagator 

gi[(m(t)-j)r+s(t)]L£ converges to the projection Pq^ = as t ^ cx), in the weak 

sense on Tig . Define the orthogonal projection = 1 — Qp on TL, where Qp = 
^^^I'djiPug). We want to show that 

lim (n, (B'YB'Un e'^^i • • • e'^^f e^^^^'+i • • • e'^^-e^'^^+i 

t—^0O 

To do so we split the Hilbert space as 7i = Hi H2, where Hi = Hs ®m>p+i '^S-, and 
H2 = ®l^iH£, and set 

0^?^+i(So))e-'^^-+ie-*^^- • • • e-'^^-+i [VLs ®m>p+i ^s] e Hi, 

= g^(®^Lie^[("^-^)^+^]^^Ai) ®?=if^f e7^2. 

Qp is a sum of terms, each one containing the operator Pq^ = Ig — Pq^ acting on at 
least one of the p factors in H2- Consequently we have 1^2 ^0, weakly in H2, as 
t ^ 00. Since ip^^^^ is uniformly bounded in t it follows that ip"^^^^ ip'2'^^^ converges 
weakly to zero in H, as t ^ 00. This proves relation (4.13). 

Thus, in the limit t — ^ 00, the only contribution to (4.12) comes from the part 
where all the free propagators e'[('"~*)'^+*l^£ are replaced by Pq^. This shows that 



lim 

t— >oo 



E{t) - {n, [B'YB'Un e'^-^i • • • e'^-^^e'^-^^'+i • • • e'^^^-e'^-^^+i 



p 



xe--Wie--^'« . . . e--^^+iO> J] U (Bj)^^ = 0. (4.14) 

i=i j=i 

We may now, as we did above, turn the operators Lj in (4.14) into K^s, also for the 
remaining indices j = p + 1, . . . ,m + 1, to arrive at 



lim 

t— >oo 



E{t) - {n, [B'YB'Un e'^-^i • • • e^^-^^Pive'^-^^+i • • • e'^-^-e'^-^^+i 



xAs i?„^+,+i(r| ' ^^"^\Bj)) ® ^m+i[Bo)n) 

p r 

i=l j=l 

where we introduce the projection Pn = 'S>m>N+i'&m[Pn£) (that projection comes from 
the left factor of the inner product and slips through [B'Y B'Un and through the first 
N propagators). 
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We have 

^m+i{Pae)e''''-+'{As(^^m+i{Bo))n = Do{s)n, (4.16) 

where Dq(s) is a linear operator acting nontrivially only on Tis- Do{s) depends on ^5, 
Bq, the interaction V and s = s{t), but it is independent of m = m{t). In the same 
way we define Di{s) G B{Hs) by 

^m{PasW^''"'^m{Ti{B.i))Do{s)n = Di{s)Do{s)n, (4.17) 

and then 1)2(5), ... , D^{s) G B{Hs)- Hence the inner product in (4.15) can be written 

as ^ 

'n, {B'YB'Un e'^^^ ■ ■ ■ e^^-^'^PTve'^-^^+i • • • e'^^^-W{s)n) , (4.18) 



where D{s) = D^{s) ■ ■ ■ Do{s) G B{Hs)- Since = P^l, where P is the projection onto 
ric, c.f. (1.8), we have 

Pive'^^^^+i • • ■ e'''^"^-W{s)^} = Pe^^^^+i • • • e'''^'^-^PD{s)n. (4.19) 

The reduced product of the propagators on the r.h.s. is the product of the reduced 
propagators, as we show in the following proposition. 

Proposition 4.1 For any q > I, let t\,...tq G and let mi, . . . ,mq> 1 he distinct 
integers. Then we have 

Pgi*l^"il . . . Q'iq^mqp _ pQ^tlKm^ p _ _ _ p^HqKmq p 20) 

Proof of Proposition If Q is a projection we set Q-^ = 1 — Q. We have 

pJ-^itqKmq p g jiqj^ (-ff^) sluce ou all factors of He with label m / mq the projection 
'!?m(-Pf ) coming from P commutes with e'*^-^'"? . It follows that Pe'*!^""! • • • p-Lei*?'^^'"? p = 
p^itiKm^ • • • ??^^(P^)P-^e'*«-'^'"9P = 0, because i?mp(-P/) can be commuted to the left 
to hit P. Thus we have Pe^*i^™i • • • e^^i^^^q p = Pe'^^iKm-, . . . Q^tq-iKmq.i p^itqKmq p ^^^^ 
we can repeat the argument. This proves the proposition. ■ 
According to (4.20) the r.h.s. of (4.19) equals M"^-^'^ D{s)n, where M = PMP 
is the operator introduced in (2.18). We then obtain the following result from (4.19), 
(4.18) and (4.15): 



lim 

t—^00 



E{t) - (n, (B'YB'Un e'^^^ ■ ■ ■ e'^-^^PM"^W-^-^PL>(s(t))0 



p r 

^II<^^)n.n(^^V|=0- (4-21) 

i=i j=i 



In order to further simplify the scalar product in (4.21) we use the ergodicity assumption 
(E). We have M^W"^"^ ^ tt = as t ^ oo, in the topology of B{Hs) (c.f. 

(2.19)), and since D{s{t)) is uniformly bounded in t we obtain 



lim \Eit) - {n, {B'yB'n) {n*s, PDisit))pns) H (a,)^^ [] {b,)^^ 

i=i j=i 



0. (4.22) 
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Remember that tp is approximated by B'Q, i.e., given an arbitrary e > we choose B' 
s.t. \\ijj-B'n\\ < e. Thus 

{n, {B')*B'n) = \\B'nf = (HV^H + 0(e))' = l + 0(e), (4.23) 

which we can use in (4.22) to arrive at 



Mm 

t—^oo 



E{t) - {n*s, PD{s{t))Pns) n {A^)ne n (BjK I = 0- (4-24) 

i=i j=i 



Note that if 1 is a degenerate eigenvalue of M, (4.22) becomes 

p 

lim 



t^oo 



E{t) - {n, {B'rB'7rPD{s{t))Pns) U (^i)n, H (^^K = «■ (4-25) 



2£ 

i=l j=l 



Thus. E{t) has a r— periodic asymptotic behaviour (s(t) is periodic) but which depends 
a priori on B' and thus on the initial state u: the simplification due to (4.23) does not 
hold anymore. 

Finally by the definition of D{s), c.f. (4.18), (4.17), (4.14) we get 

{n*s,PD{s{t))Pns) (4.26) 
= (n*s, Pe'^^1 • • • e^"^^e'^^«+i [As t|^-')^+'(B_^) ® • • • r|(5_i) Bq] P^s 



Using the invariance Pfis = e '*^<+ie '"^^^ • • • e '"^^iPfi^ and passing from if's to L's 
we obtain 

(f]J,PD(s(t))PJ75) 

= (OJ, Pe^'^i • • • e'^^'e''^'+' [As r^^"^^"+'(B_^) • • • 

• • • r|(5_i) 5o]e-^"^^+ie-^^^^ • • • e-'^^'Pn. 
= (^n*s, Pe'^^i • • • e'^-^^e'^^^+i {U^)* [As^B^^^---® Bq] 



= (n*s,PaZi^'^*HAs<^B_e(^---(^Bo)Pns) . (4.27) 



We conclude the proof of (2.29), (2.30) and hence the proof of Theorem 2.4 by plugging 
(4.27) into (4.24). ■ 

4.2 Proof of Theorem 2.5 

The proof is a simple modification of the proof of Theorem 2.4. We indicate the main 
steps. Proceeding as in Section 4.1 we see that (compare with (4.12)) 

C{t; A, As, Bo) = (n, {B'fB'A Un^^^' ■ ■ ■ ^^^rn^^'^m+ij^^ ^ ^„+i(So)o\ , (4.28) 
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where we may assume that A acts trivially on factors of TCc with index > N. As in Sec- 



tion 4.1 we replace the product of the propagators in (4.28) by e' 
Taking the limit m = m{t) ^ oo then yields 



■■e' 



lim 

t— >oo 



C{t;A,As,Bo) - {n,{B'yB'An) u;+{Pa'^\As®Bo)P) 



0, 



(4.29) 



with a speed of convergence dictated by Proposition 2.2. To complete the proof of 
Theorem 2.5 we notice that {n, {B')* B'AQ.) = {A}^,^ = u;{A) + 0{e), for arbitrary 
e > (c.f. the argument before (4.23)). ■ 



4.3 Proof of Proposition 2.1 

Let Kjn, Ljn be the operators on H that act trivially on all factors except on Hs times 
the m-th He, where they act as K, L. Given any t eM., A,B e dJls (C 9Jl) we have 

{BQs ® nc, e'*^i • • • ^tLmj^Q-itLm . . . Q-itLi^^ ^ 

= {BQs ^ ^c, e'*-^^ • • • e'^^-^Os (g) nc)-^ 

= {Bns,{Pe''''PrAns)^^, (4.30) 

where we use that K and L implement the same dynamics on 9Jl and (2.16) in the first 
step, and Proposition 4.1 in the second step. Since B^s is dense in Hs it follows from 

(4.30) that 

WiPe'^^PrAflsW < \\A\\. (4.31) 
So far we have not used that dim 7^5 < oo. In the finite-dimensional case, Tls^s is 
not only dense in Hs, but for any ^ G Hs there exists an ^ G Tls s.t. V = A^ls- Thus 

(4.31) and the uniform boundedness principle give that 

sup ||(Pe"-^P)"^|| <oo. (4.32) 

teS.,m>0 

The facts that the spectrum of Pe^^^ P lies in the unit disk in C, and that all eigen- 
values on the unit circle must be semisimple follow from the uniform boundedness of 
II (Pe'*'^P)'"|| in m. They can be shown using an easy Jordan canonical form argument. 



4.4 Proof of Proposition 2.6 

Given any normal state uj ofdJl and any unitary U onH wc have the following relation 
[7] (our definition of entropy differs from the one in [7] by a sign) 

Ent(w([/* • U)\uJo) - Ent(a;|a;o) 

= u (u* [Pe J2 + PsLs] U-PeJ^ L£,k - PsLs) ■ (4.33) 

k k 

In case J7 is a dynamics of the system, (4.33) is called the entropy production formula. 
We take U = C/Ri(t). The argument of uj in (4.33) can be written as 

f3s{aii ( J2 Le,k + Ls)-Yl -Ls}-{Ps- Ps) {a^Ls) - Ls) , 
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so it suffices to prove (2.43) for Ps = Ps = P- We want to show that 

■ ^ Le,k + Ls)-Yl - ^5 = A^(i) - akI(K^(t)+l) + ^i- (4-34) 



"RI 



It is clear that the sums in the l.h.s. of (4.34) extend only from A; = ltoA; = m + l. 
We examine the difference of the first two terms on the l.h.s., for A; = m + 1. With 
t = mr + s, we obtain for this difference the expression 

— OiRl{L£,Tn+l + Ls + V„i+i) — aRi(T4ii+l) — -^f.m+l 

= a^{{Ls)+a'^{{Vrn+i)-aUVm+i) 

= a^{{Ls + Vm) + am " K^) - ala{Vm+i) 

= a^{{Ls + Vm)+j{m)-ai,{Vm+i), (4.35) 

where wc use in the second step a^[{L£^rn+i) = -^f,mH-i) cind in the last step wc use 
definition (2.35). We now add to (4.35) the expression a^j{Lg^m) — Lg^m = ct^iiLe^m) — 
L£,m (i-e., the term with A: = m in the sums of the l.h.s. of (4.34)) and repeat the 
manipulations leading to (4.35) to obtain for this sum the expression 

a'^-^^^{Ls + y„-i) + j{m - 1) + Jim) - a^i(y„+i). 

It is now clear how to continue this process until all terms with A; = l,...,m+lin the 
sums of the l.h.s. of (4.34) are taken care of. We obtain 

m(t) 

l.h.s. Of (4.34) = ^ j{k) + Fi - aki(Hn(t)+i)- (4.36) 
k=i 

Definition (2.41) yields the result (2.43). 

To see that Ent(a; o ckj^jIcjo) is continuous in t > we show that 

Ent {u o a™f |wo) - Ent (w o a^^f '^^^^'' {loq) (4.37) 
converges to zero, for all m > 1, as s J, and s' t r. Expression (2.43) yields 

(4.37) = Lo(^P£[AE{mT + s)-AE{{m-l)T + s')] (4.38) 

-a--+^ {PeVm+i + ApLs) + ak7"'^"^'' {PsVm + A/3L5)) , 

where A/3 = Ps — Ps- The first argument of u in (4.38) equals Psj{m), see (2.41). 
Next, 

<r^(wi) - 4r'^"^''(H„) <f (vw+i - Vm), (4.39) 

as s 1 0, s' I r. But the r.h.s. of (4.39) is just j{m), see (2.35). This will compensate 
the contribution of the first argument in lu of (4.38) in the limit. Finally, a^^'^'^ (Ls) — 

•^Ri ^''^^^ (-L5) tends to zero, as s | 0, s' t r. Consequently, (4.37) vanishes in the 
limit. I 
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4.5 Proof of Proposition 2.7 

By (2.43) and (2.41), the change of entropy within an interval of duration r is given by 

Ent(a; o a|^^|u;o) - Ent(a; o a^ilwo) (4.40) 
= ^fa;(j(m + l)-a*+^(VW+2) + afii(V;n+i)) (4.41) 
-Ws - Ps)oj{a'+,^Ls) - aULs)), (4.42) 

where j{k) is given in (2.35), and m = m(t). Taking into account (2.35), (2.38) and 
Theorem 2.3 we see that the hmit t ^ oo of (4.41) is P£U}+{j+). We claim that (4.42) 
vanishes as t ^ oo. This can be seen in the following way. An application of Theorem 
2.4 shows that 

lim a;(a^V(L5) - aULs)) = T a.+ (Paf,i(i[y, L5])P)ds, (4.43) 

t-»00 Jq 

where we use that [V, Ls] eMs^Ms, which follows from the fact that e'^^^Ve-'^^^ G 
Ms (8) Ms for all t G R. Therefore 

Ent(a;oa^y|u;o)-Ent(cc;oa^i|u;o) uJ+{j+) - {Pe - Ps) /^(^"Rilii^, ^5])^)ds 

Jo 

(4.44) 

as t — GO. On the other hand, (2.46) shows that 
1 

- [Ent(a; o a^+^l^^o) - Ent(u; o a*Ri|a;o)] ^a;+(j+), (4.45) 

in the limit t ^ oo. Then the following general fact proves that (4.43) must be equal 
to zero: 

If a locally bounded function / on M has the property f{t+T) — f{t) — > a, as t — ^ oo, 
for some a G M and some r > 0, then f{t)/t a/r, as i — > oo. 

This concludes the proof of Proposition (2.7). ■ 

4.6 Proof of Theorem 3.1 

Using a Dyson expansion we get 

girx^ = e'^-^o + iA rdte^(^-*)-^°H^e"^° (4.46) 
Jo 

-A^ r [ e'^^-*^^°We'^*-'^^°We''^" ds dt + 0{X^). 
Jo Jo 
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Inserting (3.5)-(3.14)-(3.16)-(3.17) in (4.46), one gets after a somewhat lengthy but 
straightforward computation 

2 



Ma = e^^-^^ +iAsin(r) 



-l3h/2 



9P 



ff 

Jo Jo 



R+JM.+ 



(n)f ||ff/3(r2)fe-^'-ie^(*-''-)(''^-'-^)dridr2 



^0 



_gi(T-2t+2s)a. ^ e-'^^^) dsdt 
Il5/3(ri)f ||5/3(r2)fe-^'"^ei(*-^)('-^-'-i)dri dr2 

+ 

_^||g-/3h/2^^||4 



(4.47) 



Using perturbation theory [8], we then find that Mx has 4 distinct eigenvalues: 1, eo(A), 
e+(A), e_(A) which are given by 

eo(A) = 1 - AV2(ai + as) + 0{X^), (4.48) 



e+(A) 



,2ir 



AV^ 



(ai + 02) + iA^T^^IIe "■"''ff/jH^ 



//<- 



-/3/i/2^^||4 

-/3n+e-/^^^)||5Mri)||^||5Mr2) 



(4.49) 



X 



1 — sinc(r(ri — r2 — 2)) ^ 
''I - ?'2 - 2 



+ 0(A3), 



e_(A) = e' 



-2iT 



1 - 



AV^ 



(ai+a2)-iAV(||e-^'^/Vll^ 

- j j +^-^''mm{n)\\l\\9p{r2 
1 — sinc(r(ri — r2 — 2)) 



(4.50) 



, +0(A3), 

^1 - r2 - 2 /J 

where the a^-'s are defined in (3.21). Since they are strictly positive numbers, this 
proves that for |A| small enough, the operator Mx satisfies Assumption (E). 

It remains to prove that the asymptotic state 0J+,x is indeed given by (3.20). For 
that purpose, it suffices to compute the eigenvector ^^(A) of for the eigenvalue 1, 
which is obtained by perturbation theory: 



\/2 , , a2\/2 , 
^sW = - — r— V'li + - — r— V'22 



Ctl + 02 



ai + a2 



(4.51) 



A ^Pu + i^2i) + o(A2). 

V2(q;i + a-2) 
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4.7 Proof of Theorem 3.2 

Let Pq denote the spectral projection on the kernel of Lq, and fig the main term in 
the expansion of J7J(A), i.e. f^J(A) = + 0(A) (see e.g. (4.51)). Then, using (2.23), 
(2.38) and perturbation theory, one gets 



Jo 



-itLo _ JtLo 



-1 / (OS ® I WPo^^^° V^s ® Of) 

'0 



+ 0(A3). (4.52) 



Moreover, as a general fact, one has Pq^o <8) O5 = Oq (g) Og . Now, the quadratic 
interaction satisfies PqWPq = 0. Hence the second term on the right hand side of 
(4.52) cancels. 

Using (2.50), (3. 16) and (4.51) we thus have 



dS. 



+ 



I \\9p{r)f\\9(3{r')f{2-r + r')smc' (^^ 



x{a2e-^'' - aie-^''')drdr' + O(A^). 
Inserting the expression for aj in the integral, one finally has 



dS. 



+ 



2(ai + a2 



I ll5Mn)f||5Mr2)|p||5Mr3)lPll5Mr4) 



(4.53) 



xsmc 



2 /t(2- n +r2)\ . 2fTi2-r3 + r4) 



smc 



: (r2 + rg - ri - u) (^e-^(^i+^^) - g-^^^^+rs)^ + o(;s,3)^ 



where the integral runs over the four variables rj G M"*". The result follows then from 
the fact that, for any real x and y, {x — y){e~^y — e~/^^) is non-negative. ■ 



4.8 Proof of Theorem 3.3 

The proof goes exactly in the same way as for Theorems 3.1 and 3.2. We just give 
the expressions for the eigenvalues (different from 1) of Mx as well as the one for the 
entropy production. 

eo(A) = 1 - \\^{ai + as) + O(A^), 



e+(A) 



,2iT 



2^2 



AV 



(ai + a2) + iA^r^ 



JUr 



(4.54) 
(4.55) 



1 — sinc(T(2 — r)) ^ 1 — sinc(r(2 + r)) 



2-r 



2 + r 



+ 0{X% 



e+(A) = e 



-2iT 



2^2 



1- 



AV 



■{ai + 02) — iA^r^ 



/ Mr) 



(4.56) 



1 — sinc(T(2 — r)) 1 — sinc(r(2 + r)) 



2-r 



+ 



2 + r 



+ 0(A'), 
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dS-, = ^^/lb/3(r)II^II^M0ll^mc^(^^) (4.57) 

-2 (^^^^) (r + rO(l-e-/^(^+'^')) 



smc 



2{ai + a2) 



X 



2 



r(r-2)\ . 2/^-^(^'-2) 



sine smc 



2 \ 2 



+smc2 y 2 ) '''''' [ 2 )\ + 

All the integrals run over M"*", and in the formula for the entropy the integration is 
computed w.r.t. both r and r'. ■ 

4.9 Proof of Theorem 3.4 

Once again the proof goes the same way, and we only give the expressions for the 
eigenvalues (different from 1) of Mx and for the entropy production. 



e+(A) 



2^2 



1 - 



2(l+e-/3^£) 



(ai+a2 + (l + e-^-^£)sinc2 



+i 



x(|ap + |d|2-ad-ad)) 



AV^ 



+(l-e- 



smc 



tE£ ^ ad — ad 
2i 



^ t{Es + Es) ' ' 



+o(A'), 
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e_(A) = e-'"-^^ 



x(|ap + |dp 



-ad- 

sinc(r^^l^|2 _ 



ad)) 



+(1 - e-Z^^^sW (^)^^ 
r(ii/£ - Es) 

_n + „-/3E.^ l-smc(T(£;£ + £;s)) , 



+ 



A^/3r.£;g(l-e-^-^g) 
(ai + a2)(l + 6-/3-^0 



sine 



+2|bp|c|2(l + e-^^^)sinc^ j^r(£^+£5)^ 

"'''+0(A3; 



X sine 
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